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ABSTRACT 

Pinning of superfluid vortices to magnetic flux tubes in the outer core of a neutron 
star supports a velocity difference of ^ 10^ cm between the neutron superfluid 
and the proton-electron fluid as the star spins down. Under the Magnus force that 
arises on the vortex array, vortices undergo vortex creep through thermal activation or 
quantum tunneling. We examine the hydrodynamic stability of this situation. Vortex 
creep introduces two low-frequency modes, one of which is unstable above a critical 
wavenumber for any non-zero flow velocity of the neutron superfluid with respect to 
the charged fluid. For typical pinning parameters of the outer core, the superfluid flow 
is unstable over wavelengths 10 m and over timescales of ~ (A/1 m)^/^ yr down 
to ^ 1 d. The vortex lattice could degenerate into a tangle, and the superfluid flow 
would become turbulent. We suggest that superfluid turbulence could be responsible 
for the red timing noise seen in many neutron stars, and find a predicted spectrum 
that is generally consistent with observations. 
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1 INTRODUCTION 

The dynamical role of hydrodynamic instabilities in the liquid core of a neutron star is of considerable interest in understanding 
a variety of phenomena that includes spin glitches, stochastic spin variations ("timing noise") and thermal evolution, as well 
as possible precession and r-modes. One instability that could arise is the "Glaberson-Donnelly" counterflow instability 
observed in superfluid liquid helium, w herein flow of the fluid 's normal component along the vortices that thread a rotating 
superfluid drives the system turbulent IjGlaberson et al.lll974l ). Early in a neutron star's evolution, when the temperature of 
the core is comparable to the condensation temperature of the neutrons, differential rotation in the core, resulting from a 
glitch or possibly causing it , drives an Ekman fl ow al ong the rotation ax i s that can excite the Glaberson-Donnelly instability 
i Peralta et "Zlboosl . lioOQ ). IPeralta et al.l l|2006l ) and lMelatos &l Peraltal ([20071) have shown that transitions between laminar 
flow and fully-developed turbulen ce could drive spin glitches. Unstable shear layers in the outer core may also play a role in 

glitches jPeralta fc MelatosllioogI ). 

The protons of the outer core are predicted to form a type II superconductor ( Migdaj 19591 ). When the protons condensed 
early in the star's life, the expulsion time for magnet ic flux was far gre ater than the nucleation time of magnetic flux tubes, 
and the magnetic flux became confined to flux tubes I Bavm et al. 19691 ). The magnetic field inside a flux tube is comparable 
to the lower critical field for superconductivity. Be ~ 10^^ G. Entrainment of the pro ton and neutron m ass currents causes 
magnetization of the neutron vortices with a field in the vortex core comparable to Be ("Alpa r et al]|l984bl ). If the vortices are 
free to move, momentum exchange between the charged fluid and the neutron superfluid through the scattering of electrons 
with the magnetized vortex cores quickly damps velocity differences between the two fluids. The strong magnetization of the 
vortex cores, however, will pin the vortic es to the flux t ubes with a pinning energy of > 100 MeV per intersection (see [Jones 
[l991: Mendel] 1991; Chau ct al. 1992; Ruderman et al.. 1998, and the estimates of §2 below). Pinning plays a crucial role in 
the dynamics of the core fluids; to the extent that pinning is perfect, the motions of the two fluids are largely independent, 
apart from modiflcations to the average mass currents by the entrainment effect. 
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Siderv et al. I l|2008l ) and Iciampedakis et alj (|2008bl lah have studied a variant of the Glaberson-Donnelly instability for 



the core mixture of superfluid neutrons, superconducting protons, and electrons. The proton-electron liquid acts analogously 
to the normal component in liquid helium. They found that inertial waves of the neutron superfluid are subject to the 
Glaberson-Donnelly instability if there is a relati ve flow between the proton-elec tron fluid and the neutron fluid directed along 
the rotation axis, whether pinning is negligible (Sidcr v et al. 20081 ) or perfect ( Glampedakis et al .!'2008b'. Q). Such a relative 
flow could result in a precessing neutron star in which pinning of vortices to flux tubes is effective, van Hoven fc LevinI ( 20081 ) 
have shown, however, that the enormous magnetic stress in the charged fluid suppresses the instability for relative flow speeds 
less than the hydromagneti c wave speed; the instability is unlikely to occur, even in a precessing neutron star. 



Andersson et al.l (|2004l ') have shown that relative motion between the neutron fluid and the charged fluid can be sus- 



ceptible to a two-stream instability, though pinning was not considered. This instability is similar to the Kelvin-Helmholtz 
instability for relative flow of tw o fluids separated by an interface, but occurs when the fluids coexist in the same volume. 
Glampedakis fc Andersson ( 2009[ l have identifled a instability that can occur in the neutron-proton mixture of the core asso- 



ciated with inertial r-modes, assuming perfect pinning of vortices to flux tubes. The instability o ccurs if the relative angular 
veloci ty difference between the proton and neutron liquids everywhere exceeds a critical value. iGlampedakis fc Andersson 
propose that this instability triggers glitches, though it has not been explained how the the critical angular veloc- 
ity difference for instability reaches its critical value nearly everywhere in the star at once, or how a glitch is subsequently 
triggered. 

From the standpoint of building realistic rotational models of neutron stars, a crucial question is if laminar flow in the 
star is generally stable or unstable. As a flrst step towards answering this question, we studied the stability of the vortex 
lattice , pinned to the nuclei of the inner crust, in the presence of the superfluid flow sustained by pinning as the star spins 
down (|Linkll2012l '). For the situation of imperfect pinning that would arise as a result of, for example, vortex motion through 
thermal activation or quantum tunneling, two low-frequency modes appear in the system, one of which is unstable over length 
scales of less than ~ 10 m, and over timescales possibly as fast as ~ 100 d. W e suggested that the instability could lead to 
turbulence similar to grid turbulence seen in liquid helium ( Smith et al. IQQ^). In this paper, we extend the analysis to the 
neutron star core and flnd an instability that is analogous to that found in crust. As the star spins down, the Magnus force 
on the pinned vortices exerts stresses on the neutron and charged fluids that drives the flow unstable, possibly over timescales 
as short as days. 

In §2, we estimate the angular velocity difference between the neutron fluid and the charged fluid that can be sustained by 
pinning in a spinning-down neutron star. In §3, we examine the stability of this state of differential rotation, and show that it 
is unstable. In §4, we estimate the spectrum of stochastic torque that would arise if the fluid becomes turbulent, and compare 
to measured power spectra in radio pulsars and magnetars. Our chief conclusion is that the superfluid in a spinning-down 
neutron star is generally turbulent everywhere that there is signiflcant pinning, hence, throughout the outer core and inner 
crust. 



2 PINNING OF VORTICES TO FLUX TUBES 



We begin by establishing the approximate density range of type II superconductivity, where the magnetic fleld will be conflned 
to flux tubes and the pinning of vortices to flux tubes will be important. The proton fluid will be a type II superconductor 
when the proton coherence length and the London length A* satisfy 

^p<V2A,. (1) 



The London length is (|Alpar et al.lll984bl ') 



A. = 30 



-Xp Pj_4 



fm, 



(2) 



where TUp is the bare neutron mass and m* is its effective mass in the medium, x. 
of protons and neutrons, and p n is the total mass density in units of 10^* g cm 
19761 : IChamel fc HaenselboO' 



p = pp/pn, Pp and Pn are the mass dens ities 
^. In the outer core, rrip/rap ~ 1/2 ((Sioberel 



The proton coherence length is (|Mendelllll99ll ) 

ip = lexyVJf ^ Ap(MeV)-' fm. 



(3) 



where Ap is the proton pairing gap. The core begins at a density of pa ~ 1.5 jRavenhafl et al.lllisi '), and the proton gap Ap 
has a typical value of ~ 1 MeV around nuclear density pi4 = 2.8 (|Elgar0v et al.lll996l ). Combining eqs. ^ and ^ gives 



0.3 



rrip/mp 
0.5 



-3/2 



V 0.057 



5/6 



5/6 



Ap(MeV)' 



(4) 



© 2012 RAS, MNRAS 000.[TlfT2l 



Instability of Superfluid Flow 3 



The protons thus fo rm a type II superconductor in the outer core. Above nuclear density the proton gap begins to faU 
( Elgar0v et al. 19961 ). and a tran sition to a type I sup erconductor occurs at several times nuclear density. Polarization effects 
may lower Ap by a factor of 2-3 ( Schulze et al.lll996h . restricting the type II region to somewhat lower density. 

Vortices and flux tubes interact with one another and pin as a result of the interaction between the neutron and proton 
conden sates. For neutrons flowing at velocity d„ and protons at velocity Vp, the mass currents take the form (|Alpar et al 
1984bl l 



Vp + pnpV 



(5) 



The mass current of each species is generally in a direction different than the velocity of either species. This entrainment 
effect is fundamentally non-dissipative. 

For a proton mass density pp and neutron mass density p„, the coefficients for the mass currents are given by 



ppp — pp 



Pnn — Pi 



pnp — ppn — Pp 



— Pn 



5m* 



(6) 



where m* is the effective mass of the neutron; Sm* = m* — nip and (5m* = m* — m„ are the contributions to the effective 
masses due to interactions in the medium. The neutron and proton mass densities, and the total density, neglecting the 
electron mass, are 



Pp — Ppp ^^ Ppn 



Pn — Pnn Pnp p — Pn pp • 



(7) 



In t he outer core, 6m^/mp ~ —1/2 and m^/mn — 1 ( SiobertJ 19761 : Chamel fc Haensej 2006|). an d ppn is negative. 

Ruderman et al. I (| 19981 ). using Ginzburg-Landau theory and the results of lAlpar et al.l (|l984bl '). have calculated the energy 



per unit length L of a superimposed flux tube and vortex line, minus the energy for infinite separation: 

Eo TT / $0 

Ppp 



So ^ TT /_£o_y ^2 rnp_ppn^^ I A, 
L ~ 8 V TT Al ) * m„ Ppp \ ^ 



= ^B^B^Al In ( ^ 



(8) 



where $o = hc/2e is the fiux quantum, _B$ = $o/7i"Aj is the characteristic magnetic field in a fiux tube core. By = 
{^Q/i^A?i,){mpPpn/mnppp) is the field in the vortex core, and ^„ is the neutron coherence length. Since ppn is negative, 
E/L is negative if the vorticity is parallel to Bv, and positive if anti-parallel. The interaction is primarily magnetic. Both the 
flux tube and the vortex line have their magnetic fields screened over the London length A, . 

If the angle between the vortex line and the flux tube is 9, the overlap length is / ~ 2A, /sin^. The interaction energy 
per vortex-flux tube junction is approximately 



Ep{e) ~ I 



Eo 
L 



B* (A,/) In 



4 \tvAV 



A: 



o (5m* 
^ ^ In 



cote 



(9) 



This equation provides only a rough estimate for arbitrary 9, since eq. ((H)) does not account for modiflcations of vortex and 
flux tube structure for 6 7^ 0. The pinning energy is approximate ly the magne t ic energy density in the overlap r egion, times 
the overlap volume A^l. Similar estimates have been obt ained by Jones (Il991 ). Mendelll (1991), and Chau et al.l (1992). 

The magnetic fleld is expected to be highly tangled ( Ruderman et all 1998 ). and so 9 could take a broad range of values. 
For an average radius of curvature of fiux tubes of the maximum overlap length is 

2A. 



= 2i?$ cos"^(l - A,/Ri,) = 



(10) 



where 6'niin is the minimum angle at which a curved flux tube intersects a straight vortex. For i?$ >> A,, this equation gives 
Smin — (2A,/ii#)^''^. The interaction is repulsive for 7r/2 < 9 < tv, but the magnitude of the pinning interaction is the same 
for both 9 and n — 9. Vortices will be immobilized by fiux tubes in much the same way regardless of the sign of the interaction. 
We average \Ep(9)\ over 9 to obtain the average pinning energy for i?# >> A, 



(Ep) 



(tT - 26»,„in) 



d9\Ep{9)\ 



■ Ep{cot( 



1) In 



2X7 



(11) 



Unless Rg, is much larger than ~ 10"^ A,, the logarithmic factor is not significantly different than one, so we henceforth set this 
factor to unity. If R<s> is much larger, however, it is probably a better approximation to set 9 equal to the inclination angle 
between the star's magnetic moment and the rotation axis. In either case eq. (|lip. with the logarithmic factor set to unity, 
should be a good estimate. 

The pinning energy per vortex-fiux tube junction is, taking a typical value ln(At/^„) — 0.5, 



(Ep) ~ 10' 



m*/mp 
0.5 



-1/2 



I (5m* I / r, 

a5~ 



(12) 
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The pinning force is Fp ~ (Bp) /A,, typically ~ 1 MeV fm^^. 

For a single vortex immersed in a tangle of flux tubes, the average length between intersections will equal the average 
distance between flux tubes, lg> — n^^''^, where n$ = -B/$o is the areal density of flux tubes and B is the average magnetic 
field strength. The pinning force per unit length is fp — Fp/Li. For a velocity difference w between the neutron superfluid 
and the pinned vortex, the fluid flow exerts a Magnus force per unit length pnnw, where k = h/2mn is the vorticity quantum. 
The critical velocity difference Wc that can be sustained is given by 

fp = -i^r" = P^i^'^c (13) 
SO that 

(S) (^) 

where B12 is the strength of the field in units of 10^^ G. The critical velocity is not a strong function of density. 

As a star spins down, pinning prevents the neutron superfiuid from corotating with the charged fiuid, and a velocity 
difference w builds. If w ~ at some time, spin down of the crust will cause the critical velocity to be reached in a time At 
that is short compared to the star's age: 

where 1/ is the spin frequency of the crust, t^gc = v j^ii is the spin-down age, and R is the stellar radius. Hence, the angular 
velocity difference between the superfluid and the charged components will be Wc/R ~ 10~^ rad s~^ throughout the star's 
life (except, for example, right after a glitch). We now examine the stability of this state of differential rotation. 



3 STABILITY ANALYSIS 

The problem of the coupled dynamics of the neutron and proton fluids can be treated with hydrodynamics for length scales that 
are large compared to the inter- vortex spacing restricting the treatment to wave numbers kl-u << 1. In a single-component 
superfluid, the vortex lattice supports Tkachenko modes of speed ct = (fi.Q./'im) (|Tkachenkolll966al lbl). where Q, is the spin 
rate of the superfluid and m is the mass of the fundamental boson, twice the neutron mass fo r a neutron superfluid. (For 
stu dies of Tkache nko modes in neutron stars and the effects of dissipation on these modes, see Noronha fc Sedrakian 20081 
and Haskell 2011 ). For a typical neutron star, ct — 10~^ (fl/lOO rad s~^)^^^ cm s~^. The degrees of freedom of the vortex 
lattice can be ignored to a very good approximation when the speed of the background flow with respect to pinned vortices is 
much larger than ct (Link 20121 ). as is the case for vortex pinning to flux tubes. The areal density of vortices is = 2mO,/h 
for a uniform vortex lattice; hence, the requirement that klv << 1 is equivalent to CTk « Q. 

In the labor atory frame, the acceleration equations for neutrons flowing at velocity Vn and protons plus electrons flowing at 
veloc ity Vp, are (|Prbj|2004l : Undersson fc Comerlbood : | iGlampedakis et al ] |2008bl : Ivan Hoven fc LevinlboO jjciampedakis et al 
2OI1I ) 



{dt + v„ ■ V)(ii„ - e„Wr, 



ap^Vi = -V/i„ + f/p„ 



(dt + Vp ■ \7)(vp + epWnp) + epw\pVvl = -V Pp - f / Pp + Ve V^Vp + (l/inpp) B^a ■ VBc 



(16) 
(17) 



where Wnp = Vn — Vp, i is a coordinate index and is summed, pn and pp are the neutron and proton chemical potentials, p„ 
and Pp are the neutron and proton mass densities (neglecting the electron mass), f /pn is the force per unit volume exerted on 
the neutron fluid by the proton-electron fluid, is the electron kinematic shear viscosity, and BcB = \/BBc, where B is the 
average field strength and Be — 10^^ G is the value of the lower critical field. Electron-electron scattering gives the dominant 
contribution to the viscosity, as all other scattering processes are strongly suppressed when both the neutrons and protons 
are superfluid. The terms proportional to e„ and ep account for entrainment between neutrons and protons, and are related 
by p„e„ = Ppep. We take the entrainment coefficients to be constants. The small effects of electron inertia and the London 



current, both neglected here, have been considered bv IGlampedakis et al.l l|201ll ). 

For the highly-conductive proton-electron fluid, the magnetic fleld is frozen to the proton-electron fluid. The field obeys 
the induction equation 



dtB + \7 X (B xvp)^0. 



(18) 



The force per unit volume on the neutron fluid, f /pn, is equal to the Magnus force per unit volume on the vortex array, 
that is, 



(19) 
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where Vy is the local vortex velocity and = \7 x v„ is the local vorticity of the neutron superfluid. For perfect pinning 
of vortices against flux tubes, = Vp. The vortex velocity with respect to the proton-electron fluid has a component along 
Wnv = Vn — Vj,, and a component orthogonal to both Wnv and The local vortex velocity can thus be written as 



Vv — Vp + (1 



jp, and a component orthogonal to both Wnp 

/?') w„p + j3 



(20) 



where /3' and /3 are coefflcients that determine the vortex mobility, and we assume to be constant. Perfect pinning of vortices 
to the charged fluid corresponds to P' — 1 and 13 = 0. The force, from eq. (|19p . is 



//Pn = /3't<J„ X Wt^p + /?a>„ X x w^p). 



(21) 



This force is the mutual friction force introduced in by iHall fc VinenI (|1956| ). but appearing here in a different context. 
Imperfect pinning, that is, "vortex creep", corresponds toa = l — /3'<<1 and /3 << 1. We refer to a and j3 as the "pinning 
coefficients". Perfect pinning corresponds to the limit a = /3 = 0, while no pinning (/ = 0) corresponds to a = 1 and 13 = 0. 
Vortices move with a component along Wnp, so that < q < 1. The energy dissipation rate per unit volume is determined 
by P, which must be posi tive to give local entropy production. More general dissipative forces have been considered by 
Andersson fc Comeil (|2006l l. 



Vortex creep could be a low-drag process, with /? << a, or a high-drag process, with /? >> q. In much previous work on 
pinning, the high-drag limit has been implicitly assumed through the following relationship between /? and /?': 

7^2 



13' 



7^/3, 



(22) 



l + 7^2 

where 7?. is a dimensionless drag coefficient. In this drag description, imperfect pinning corresponds to 7?. >> 1 so that eq. 
(|22|l requires P » a. Eq. ([22} is not true in general, since the presence of non-dissipative forces between vortices and the 
charged fluid can give P « a, a. regime of low drag (,Link.2009iy As we show below, it is the low-drag regime that is likely to 
be realized, with vortex creep being unstable in this regime. A crucial feature of our analysis is that we do not impose eq. H22p 
for imperfectly-pinned vortices. We will assume in §4 that eq. (|22|l holds only for the tiny fraction of vortex length (~ 10~®) 
that is unpinned at any instant. 

To examine the stability of the flow with imperfect pinning, we use a local plane wave analysis in the frame rotating 
with the unperturbed charged fluid at angular velocity CIq, in which Vp = and the unperturbed flow velocity arising from 
spin down of the charged fluid is wo = lUnp.o = wqx. Restricting the analysis to the regime kR^ 1, where R is the stellar 
radius, the background flow can be treated as approximately uniform. Let the rotation axis z, k, and wo all be coplanar, 
with an angle 9 between k and the rotation axis. The angle between z and B is 6b. We restrict the analysis to the quadrant 
< ^ < 7r/2. The linearized equations of motion in the rotating frame are: 



{dt + Wo ■ V){5v„ - e„Sw,^p) 
V ■ Sv„ = 



e„woV(i ■ Sv„) - 2flo X 5ii„ = -VS^i' + Sf/p„ 



dt{vp + epSwnp 



(1 



epWo\/{x 
a) S{u:„ X Wnp 



SVn) + 2flo X Vp = —VSfl'p — Sf/pp + UeV 
} + P (5{a)„ X {lJ„ X Wnp)}. 



Sf/pn 
V • = 

dtSB + V X (Bo X Vp) 

where S denotes a perturbed quantity, and 5vp = Vp. Here Pn^p = Pn,p — Pn,p(fio x r)'^ /2. The vorticity appearing in this 
equation is the total vorticity evaluated in the laboratory frame, and is only slightly larger than 2i}o'- 







(23) 
(24) 

(25) 

(26) 
(27) 
(28) 



Wo 



10-3 f ^!^) ( 

V105 cm s-V V 







100 rad s' 



(29) 



2^0 -R!^o 

We henceforth take Wo = 2f2o; we have conflrmed that the results presented below are essentially unchanged for |[Jo— 2f2o! << 1. 

For the shear perturbations we are considering, k ■ 5vn = and k ■ Vp = 0, that is, the velocity perturbations in the 
directions y and e = — cos6'i: + sin6 z are orthogonal to k. We Fourier transform (oc g^k r-icrt-j [23]-[28] and project 



the acceleration equations p3|l and (|25|l onto y and e. Defining a' = a - 
cb = cos{6 — 0b), and Xp = pp/pn, we obtain the system of equations: 

A — icr' — ienkwQS — 2cf2oQ —A + ienkwoS 

c(2f2oQ — ikwo sP) {en - 1) io-' + 2r2oc^/3 - ikwosP' C 

~2Xp^nocP' 



kwosinO, c = cos 9, s = sin 6', 



BoBc/4npp, 



B 



x-\-A + 2^os'P) 



2^00 



-~2QocP' 
'A + ienkwos + 20,08^ P 
2VLoc{x-'^P' - 1) 
B - 2x-^Q.os^p 



y ■ 5Vn 
e ■ SVn 

y-Vp 
e ■ Vp 



(30) 



where 
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A = en icy + 2fio/3 — ikwosP' , (31) 

B = —ia + ivgCgk'^ /a + Vek^ + x^^ A, (32) 

C = c{ikwosl3 + 2Q.ol3'). (33) 
(Recall that /?' = 1 — q). The speed ol hydromagnetic waves is 

1/2 / ^ ^ -1/2 



VB=2X 10** 



significantly greater than the speed too of the background flow (see eg. I14[). 

The dispersion relation that follows from eq. (|30p is an impressively lengthy sixth-order polynomial in a in which all 
coefficients from sixth order to zeroth order are non-zero; it is easies t to explore various limits by working with eq. (1301) 
directly. We first confirm a previous result of van Hoven fc LevinI 1 20081 ) for perfect pinning (a = /3 = 0), no background flow 



(wo = 0), and no entrainment (e„ =0). In this case, the polynomial factors into the form 

cr^ (cr* + Ca^ -f {-A" + 2B + AC) (T^ + BCa + B'^) = 0. (35) 
If we define 

f{a) = a^ +Aa + B ^0 (36) 
and form the combination 

a\f{af + {C-2A)f{a)a)^0, (37) 



we obtain eq. ((35|), so two of the non-zero modes are given by the the simpler quadratic expression. Reading A and B from 
eq. (|30|) gives the dispersion relation 



m 2,2 



Uflfc = 0, (38) 



as found by van Hoven fc Levin ( 20081 ) rl For — 0, the system has only damped hydromagnetic waves that travel at speed 

Vb- 

We will not present here an analysis of the full mode structure of the system, but focus on two low-frequency modes 
that appear for imperfect pinning. For small a, /?, and e„, the two zero-frequency modes of eq. p5|l become small, and we 
can obtain these modes by working to second order in a. Since wo is much smaller than vb (eqs. 1141 and [34|l. we can further 
simplify the problem by taking the limit u_b — > oo for finite uiq. We take this limit by keeping only terms that multiply v%, 
the highest order at which vb appears. With these approximations, eq. pO|l gives 



(1 - e„)^cr^ + (2elkwos + e„{-2i,3(l + c^) - 2/?'fcmos} - 2akwos + 2i/?(l + c^))a 

- {e„kwosf + e„{2p'{kwosf + 2ipnokwos{l + c^)) - 4(a^ + /3^)t^oc^ + {akwosf - 2ial3kwos = 0. (39) 

The electron viscosity does not appear at this level of approximation; in the ws — > oo limit, the flux tube array is infinitely 
rigid, and vortex motion proceeds without producing shear in the proton-electron fluid. The validity of this approximation is 
conflrmed below. The solutions to eq. (|39|l are 



a± — - — ^ — ({a — e„)k'wo sinO — iQo{l + cos^ 6)P ± {A^I^o? cos^ 9 — SIq/^^ sin* 6 — 2iaPQ,okwQ cos^ ^sin^) \ . (40) 
For 9 = Q, the modes are 

a± = -^(±a-i/3). (41) 

1 — tn 

Imperfect pinning has introduced two low-frequency modes to the system that are associated with slow vortex motion under 
the Magnus force. The modes are underdamped for /3 < a, which deflnes the regime of low-drag creep that we will explore 
further. 

The general solution with eigenvalue a- in eq. (|40|) is unstable above a critical wavenumber kc- To find this wavenumber, 
we write the solution as 



a- = A + iB - VC + W, (42) 
where 

^ Our definition of a differs from that of Ivan Hoven fc Le"vinl ll2008h by a minus sign. 
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B = 



1 — En 1 — En 

The solution becomes unstable when 
Im(a-) = B - Im (VC + iO) = 0, 
where 

Im [VC + iD) = {C + D f''^ sin ( i cos 



C = 



(1 



I 2 2 o2 4n 

-(q c -P s) 



D = 



kwoc s 



C 



VC2 + 



Combining this equation with eq. (|44p gives 

Solving for k gives the critical wavenumber fee above which the system is unstable: 



flo {0^ ^a^fl'^ l + cos^e 



(43) 



(44) 



(45) 



(46) 



(47) 



Wo Q sm y cos f 

independent of entrainment. The critical wavenumber fcc is minimized for Q = tan^^(\/2). For k » kc, we have the approxi- 
mate solutions 
1 



-((a — e„)kwo sin =p i(Q/3 Slofcwo cos 6 sin S) 



l/2x 



k >> fee 



(48) 



For the outer core, e„ ~ 3 x 10~^ for Xp = 0.05 i Chamel fc Haensel 20061 ) . Though entrainment is essential in producing 
vortex pinning, it has a negligible effect on the growth rate of the instability, and no effect on the critical wavenumber, so we 
ignore it in the estimates below. 

The instability arises from coupling between the velocity difference Wnp and the neutron vorticity a;„ through the force 
of eq. [26]. Dissipation damps perturbations for k < k^, but for k > kc the finite vortex mobility gives rise to growing 
perturbations under the Magnus force. For k >> kc, the growth rate scales as (ap-wo)^^^ . For /3 << a, kc takes a constant 
value, but the growth rate of the mode becomes small, going to zero as j3 goes to zero. In the highly-damped regime, /3 >> a, 
damping restricts the unstable mode to large k, generally stabilizing the system. There are no unstable modes for either a = 
or P — 0; the instability occurs only if vortex motion has components along both Wnp and (d)n x Wnp) x ijj„. 

We now confirm that the viscous stress is negligible for the a- unstable mode. The magnetic stress force per unit volume 
on the charged fiuid from ea. (|25p is 

(49) 



47rpp 

while the viscous force is 

fv/pp = J^eV^Vp 

Fourier transforming, and using the induction equation H28[) . gives 



uc (Bey- „ ,. 



This ratio is largest for high wavenumber. In this limit 

5fm 



(50) 



(51) 



(52) 



In the outer core, Vc is typically ~ 10® cm'^ s ^ ( Cutler fc LindblomlllQSTI : lAndersson et al.ll2005f ). Below we estimate a/3 ~ 
■j^g-18 rpjjg hydrodynamic treatment requires fecr << ^o, so &fv/5fm ^ 10~® for a typical neutron star with ilo = 100 rad s~^ 
and Bo = 10^^ G; the viscous force is negligible compared to the magnetic force for these low-frequency modes. 



4 ESTIMATES OF THE INSTABILITY GROWTH RATE 

To obtain the growth rate of the instability, we now estimate the pinning parameters a and /3 for the vortex creep process. To 
make these estimates, we regard the process of vortex creep as consisting of two distinct states of motion for a given vortex 
segment. Most of the time, the vortex segment is pinned. A small fraction of the time, the vortex segment is translating 
against a drag force to a new pinning configuration. The mutual friction force we are using (eq. I2ip is 

f / pn = aj„ X Wnp — aU)n X Wnp + (3 UJn X {u)n X Wnp)- (53) 
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This force represents the spatially and temporally averaged force exerted on the neutron fluid by the creep process. The first 
term is the Magnus force for perfect pinning, while the remaining terms give the contribution to the force due to vortex 
motion. When a vortex segment is unpinned and moving against drag, we take the force to have the same form, but with 
different coefficients: 

fo/Pn =ljJn X Wnp - Qo t^n X Wnp + A) X {lJ^ X W np) . (54) 

An unpinned vortex segment remains unpinned for a time to ~ d/wo, where d is the d i stanc e the segment moves before 
repinning. This distance is comparable to the distance between pinning sites ( Link et al. 1993l V For pinning to flux tubes, 
d ~ Z$ ~ $0 /B, and the average time that a vortex segment is unpinned is 



to ~ 10"'* B 



V 10° cm s / 

much shorter than the hydrodynamic timescales of interest. Suppose that at any instant, the volume- averaged fraction of 
vortex length that is unpinned is << 1. We now average /o over a volume that contains many vortices, and over a time 
long compared to to but short compared to hydrodynamic timescales, to obtain 

ifo/Pn.) = X Wnp — fvOLQ bJn X Wnp + fvPa iOn X {uj„ X Wnp)- (56) 

Quantities related to the flow are unchanged by the averaging procedure since the superfluid flow velocity, the vorticity, and 
the relative flow velocity are independent of whether vortices are pinned or not. The factors of fv in eq. (|56[l account for 
the f act that only the motion of the translating vortex segments contributes to the mutual friction (see, also, Jahan-Miril 
20061 ) . For vortex motion by thermal activation, /„ o-^f^ 



temperature. The value of /„ is unimportant for the following estimates. 

The force of eq. (|53p . which is appropriate for vortex creep, must equal the average force (fo/pn), giving the following 
relationships: 

a = Uao and /3 = /„ /3o ^ ^ = ^ (57) 

a Qfo 

We now take estimates of Po/ao to obtain the ratio P/a. 

The dominant drag process on unpinned vortex se gments identifled so far is the scattering of electrons against the 
magnetized cores of neutron vortices ( Alpar et al. 1984lJ ). For unpinned vortex segments, vortex motion can be described in 



terms of a single drag coefficient, as in eq. (|22|) . so that 

l-ao = Y^:^^ ^nPo, (58) 



Using recent values of the proton effective mass, Siderv fc Alpar ( 20091 ) estimate 



2 X 10"^ (59) 
giving Qo — 1 and /3o — Ti, so that 

^ = ^^7^. (60) 

Qo Q 

Vortex creep in the outer core is thus well into the low-drag regime. We stress that, in estimating Po/ao, we have assumed that 
the drag relationship given by eq. (|22p holds only for the small fraction of unpinned vortex length; the volume-averaged 
quantities /3 and a satisfy /3/a = /3o/qo, from eq. (|57p . but there is no value of TZ that gives both /3 << a and a « 1, the 
creep regime of low dissipation in which we are interested. 

We now estimate /3. We adopt polar coordinates {r,(j),z), with the unperturbed vorticity along z and the unperturbed 
ffow Wo along (j>, and take the unperturbed ffow and vortex velocity field to be axisymmetric. In the rotating frame {vp — 0), 
the unperturbed vortex velocity from eq. (|20p is 

Vv,o = ctwocp + Pwof = hv^jfi (61) 

where h is the average direction of vortex motion. 

For steady spin down of the star, the core superfiuid, the charged components, and the crust are spinning down at the 
same rate fo r a local differential velocity wo. The creep velocity in this steady state is related to the spin-down rate by 



( Alpar et a l. 1984a; Link et al. 1993,) 



n Q ■ 

n^-2—VvO-f = ~2—woP^no, (62 
r r 

where Q, is the spin rate of the superfluid, Clo is the observed spin down rate of the crust, and r is approximately the stellar 
radius R. We arrive at the estimate 
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R -,^-11 f "'o 



4it)oiagc Vl05 cms^V V 10^ yr 



(63) 



where 57 ~ f2o is assumed. Eq. (|63[) . in combination with eq. (|57p . gives the fiducial values a = 10~^ and Qf/3 = 10~^*. For 
these values, we deduce fv ~ (a/3/ao/3o)^''^ ~ 10~^, that is, most of the vortex length is pinned at any instant. 

We can now proceed with estimates of the instability length scale and growth rate. For a >> /? and 6 = tan~^(-s/2) in 
eq. (|47[) . the critical wavenumber is 

fc. ^ 6 = 6 X lO-^* f ^) ( , " , (64) 

Wo VlOO rad s^V VlQS cm s" V 

corresponding to a wavelength A = 27r/fc ~ 10 m. For k >> kc, the growth rate from eq. (|48|l is 

^M-)-2(^j il^^TIII^J 1 105 cms-0 ll^J '•''^ 
The hydrodynamic treatment is restricted to crk « Q. To estimate how high the growth rate could be, we consider a 
maximum wavenumber defined by crfcmax = 0.1 f2, where ct — 10~^ {Q/lOO rad s~^)^^^ cm s~^. The growth rate at this 
wavenumber, from eq. (|48p . is 



±M._,.„.,,.o..(^)"'(^)""(,^)- 

For f2 = 100 rad s~^, the corresponding wavenumber is fe ~ 100 cm~^. Eq. H66p does not represent a physical limit, but only 
the restrictions of the hydrodynamic treatment; the instability could continue to exist also for wavenumbers in the regime 
kcT > 57. If vortex creep is in the strongly-damped regime /3 » a, contrary to the estimates here, there is still a broad 
window for instability. Requiring kc < fcmax gives 

^ < 2 X 10* f f ^) a, (67) 

V105 cm s-V VlOO rad s-V ^ ' 

and the star will be unstable at some wavenumber that is consistent with the hydrodynamic approximation crk 



5 TIMING NOISE DRIVEN BY SUPERFLUID TURBULENCE? 



The superfluid flow could quickly become turbulent through the instability identifled here. Turbulent flow would exert random 
torques on the crust, possible contributing to observed timing anomalies such a s timing nois e and glitch es. We now estim ate 
the spectrum of turbulent torque fluctuations using the similarity argument of iKolmogorovl (|l94ll ) and lObukhovl (|l94ll '). In 
this argument, the properties of the turbulent spectrum follow from the assumption that turbulent energy is being dissipated 
at a rate e per unit mass at some high wavenumber. 

The torque fluctuation per unit mass of fluid for a mode of wavenumber k is 



SN{k) dk ~ ±1 {Av)f 



(68) 



where Z ~ 1/fc is the characteristic dimension of a turbulent cell, and Aw is the characteristic change in flow speed over a time 
1//, and fl ~ Av. The units of 3N{k) are cm^ s^^, and the only possible combination of e and k with the correct units is 



SN{k) 



(69) 



where A is a dimensionless constant. The time average of 5N{k) is zero. 

Integrated over the star of radius R, the contributions to the total torque by different cells of wavenumber k add 
incoherently, and the fluctuation in the total torque at wavenumber k scales as 

(70) 



SNtot{k) oc VM6N{k), 
where M is the characteristic number of turbulent cells of wavenumber k, 



M = 



(?) 



(71) 



Here d is the dimensionality of the turbulence. Isotropic turbulence corresponds to d = 3, whil e highly-polarized turb ulence 
corresponds to d = 1. The turbulence should be polarized to some extent by rotation (see, e.g.. lAndersson et al.ll2007l '). 
The fluctuations in the total torque scale as 



and the spectral power of torque fluctuations (spectral index at) is 



(72) 
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SOURCE SPECTRAL INDEX REFERENCE 



PSR 0823+26 


at = 


-1.16 ± 0.48 


Bavkal et al. 


(1999) 




at = 


-0.39 + 0.36 






PSR 2021+51 


at = 


-1.95 + 0.74 


Bavkal et al. 


(1999) 




at = 


-1.64 + 1.40 






PSR 1706-16 


at = 


-2.41 + 0.79 


Bavkal et al. 


(1999) 




at 


-1.94+ 1.33 






PSR 1749-28 


■ = 


-0.88 + 0.05 


Bavkal et al. 


(1999) 




at = 


-1.05 + 0.83 






PSR B1509-58 


ap = 


-4.6 + 1.0 


Livingstone et al. (2005) 


PSR B0628-28 


ap ~ 


-1 (t < 2 yr) 


Hobbs ct 


2010) 


PSR B0628-28 


ap ~ 


-5 (t > 2 yr) 


Hobbs et al. 


2010) 


SGR 1900+14 


at = 


-3.7 + 0.6 


Woods et al. 


(2002) 


SCR 1806-20 


at = 


-3.6 + 0.7 


Woods ct al. 


(2002) 



Table 1. Estimates of spectral indices. The s econd entry for PS Rs 0823+26, 2021+51, 1706-16, and 1759-29 follows from a fit to the 
arrival times to a higher-order polynomial; see lBavkal et al1ll999l for details. PSR B0628-28 shows ap ~ —1 over timcscalcs of less than 
two years, and Op ~ —5 over longer timescales. 



'Vi/3+d 



r 



to r 



d\5N,,,{f)\' 
df 

Let the rotation phase of a pulsar be (pit), so that d(j){t)/dt is the spin rate. A torque variation SNt, 
in the rotational phase 5(j>{t), where 



(73) 

t(t) produces a variation 
(74) 



and /c is the moment of inertia of the crust plus any other components of the star tightly coupled to the crust. The spectrum 
of phase residuals &()) (spectral index ctp) is therefore 

dmf)\^ 



df 



-2b/3+d 



r-' to /- 



(75) 



"Torque noise" corresponds to ap = —5. We conclude that turbulence can produce such a power spectrum only if the 
turbulence is nearly isotropic (d = 3). From the scaling argument given here, turbulence cannot p roduce "frequency noise" 
{ctp = —3) nor "phase noise" (op — —1). In early work on 45 pulsars, iD'Alessandro et al.l ([1995|) reported power spectra 
consistent with a range —1 < oip < —5. More recent estimates are given in Table 1 for six radio pulsars and two magnetars. 
These results should be regarded as illustrative only; the power spectra of m any pulsars do not follow simple power laws 
i Cordes fc Downslll985l ). and some sources show quasi-oscillatory behavior; see iHobbs et al.l l|201ol ) for a discussion of these 
issues. Though timing noise is a complicated process that is difficult to characterize, we note that the basic power spectra 
predicted by turbulent torque fluctuations are generally consistent with observed power spectra. 



6 DISCUSSION AND CONCLUSIONS 



We have identified a dissipation-driven instability that could operate in the fluid of the neutron star core over length scales 
shorter than ~ 10 m and over timescales as fast as days. The instability is a result of the forcing the vortices against the 
rigid array of fiux tubes. The superfiuid flow of the in the outer c ore could be come turbulent if dissipative processes do not 
suppress the instability. The excitation of Kelvin waves vortices (|Linkl |2004| ) will lead to dissipation, but is probably not 
nearly strong enough to put vortex creep into the highly-dissipative regime that would be required to quench the instability; 
this process damps rotational perturbations over many spin periods, and so is in the low-drag regime. Other dissipation 
mechanisms of possible relevance require study. If the core fluid beco mes turbulent, the closest experimental analogue might 
be grid tur bulence that has been well studied in superfluid helium ( Smith et al. 19931 ). In our previous work on the crust 



l Linkll2012l ). we found that the flow is unstable there as well. We conclude that superfluid flow throughout much of a neutron 
star is forced unstable by the pinning stresses on the fluid that arise as the star spins down. 

To illustrate the basic instability, we have taken a simplified mutual friction force in which the pinning coefficients a 
and P are constants. The estimates of the instability growth rate given here should be regarded as rough esti mates. For 
thermally-activated vor tex creep, the pinning coefficients will have exponential dependence on the scalar Wnp flAlpar et al 



1984al ; iLink et al.l[l993l l. We expect that this strong velocity dependence will significantly enhance the growth rate of the 



instability, without significantly changing the critical wavenumber so long as vortex motion is in the regime of low drag. 
Further work is needed on this issue. 
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Turbulent flow could create torque fluctuations with spectra that are consistent with those measured in many radio 
pulsars and in magnetars. While the simple scaling argument given in §5 suggests that turbulent torques cannot produce 
"phase noise" or "frequency noise" , it would be interesting to further explore the possibility that steeper power spectra result 
from turbulent torques. In future work we will determine if the instability identified in this paper leads to fully-developed 
turbulence, and estimate the magnitude of the fluctuations. 
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